The new form of the equation of state for dark energy fluid and accelerating 

universe 



O 

o 

(N 

(N 
> 

o 

\o 
o 
\o 
o 



Oh! 



X 



Shin'ichi Nojiri and Sergei D. Odintsovt2.|i| 

^Department of Physics, Nagoya University, Nagoya 464-8602. Japan 
''^ Institucid Catalana de Recerca i Estudis Avangats (ICREA) and Institut de Ciencies de I'Espai (lEEC-CSIC), 
Campus UAB, Facultat de Ciencies, Torre C5-Par-2a pi, E-08193 Bellaterra (Barcelona), Spain 

We suggest to generalize the dark energy equation of state (EoS) by introduction the relaxation 
equation for pressure which is equivalent to consideration of the inhomogeneous EoS cosmic fluid 
which often appears as the effective model from strings/brane- worlds. As another, more wide gen- 
eralization we discuss the inhomogeneous EoS which contains derivatives of pressure. For several 
explicit examples motivated by the analogy with classical mechanics the accelerating FRW cos- 
mology is constructed. It turns out to be the asymptotically de Sitter or oscillating universe with 
possible transition from deceleration to acceleration phase. The coupling of dark energy with mat- 
ter in accelerating FRW universe is considered, it is shown to be consistent with constrained (or 
inhomogeneous) EoS. 

PACS numbers: ll.25.-w, 95.36.+X, 98.80.-k 



I. INTRODUCTION 

The number of attempts is aimed to the resolu- 
tion of dark energy problem (for recent review, see 
Q IE 111) which is considered as the most fundamen- 
tal one in modern cosmology. Among the different 
descriptions of late-time universe the easiest one is 
phenomenological approach where it is assumed that 
universe is filled with mysterious cosmic fluid of some 
sort. One can mention imperfect fluids |^], gen- 
eral equation of state (EoS) fluid where pressure is 
some (power law) function of energy-density |^ , flu- 
ids with inhomogeneous equation of state [^j , where 
EoS with time-dependent bulk viscosity is the par- 
ticular case 0, coupled fluids [1,0, etc. The EoS 
fluid description may be even equivalent to modified 
gravity approach as is shown in p^ . 

As it has been recently discussed in it is not 
easy to construct the dark energy model which de- 
scribes the universe acceleration and on the same 
time keep untouched the radiation/matter domi- 
nated epochs with subsequent transition from de- 
celeration to acceleration. In order to minimize 
the dark energy effect at intermediate epoch one 
may speculate about sudden appearence of dark en- 
ergy around the deceleration-acceleration transition 
point. In other words, one may suppose that EoS 
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of DE fiuid is of the form p — 9{t — td)wp where 
td is transition time and w is DE EoS parameter. 
Before transition point, DE plays a role of usual 
dust which changes EoS by some unknown scenario. 
In the similar way, one can generalize other cosmic 
fluids with more complicated EoS. This introduces 
the idea of structure/form changing EoS in different 
epochs. The simplest example of such cosmic fluid is 
oscillating dark energy (12,, JL3J. Finally, the reason 
why cosmic fluid still escapes of direct observations 
could be that it has completely unexpected proper- 
ties, for instance, in EoS picture. 

In the present letter the new form of dark energy 
EoS is considered. As the first step, we introduce 
the relaxation equation for pressure (the analog of 
energy conservation law for energy-density). It is 
then shown that such constrained EoS is equivalent 
to usual but inhomogeneous EoS which is known to 
be the effective description for brane- worlds or modi- 
fied gravity 6] . The generalized inhomogeneous EoS 
which contains time derivatives of pressure is intro- 
duced. The number of examples for such EoS cos- 
mic fluids is presented and the corresponding FRW 
cosmologies are described. It is shown that cosmic 
speed-up in the examples under consideration corre- 
sponds to the asymptotically de Sitter or the oscillat- 
ing universe where accelerating/decelerating epochs 
repeat with possibility to cross the phantom barrier 
or to make transition from deceleration to acceler- 
ation. In all cases, dark energy EoS parameter is 
close to —1, being within the observational bounds 
to it. It is demonstrated that the inclusion of matter 
may be consistent with constrained EoS. 
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II. CONSTRAINED EQUATION OF STATE 

Let us discuss the possible modification of tlic 
equation of state in such a way that it would change 
its structure/form during the universe evolution. 
Consider the balance equation for the energy (en- 
ergy conservation law) 

p + 'dH{p + P) = Q. (1) 

It can be represented as a relaxation equation 

'i' = --(4'-*o), (2) 
r 

where ^' coincides with p, relaxation time r is t = 
and stationary (or equilibrium) value of p is po = 
^0 = —P- To formulate consistently this equation 
we need, as usual, the equation of state (EoS). The 
standard barotropic EoS is P = P{p), providing the 
equation for p only: 
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p=-ip + P{p)). 



(3) 



Let us conjecture now that cosmological fluid is 
described by different EoS at different epochs. In 
other words, to describe the transition from one 
epoch in cosmological evolution to another we try to 
introduce the transition from one EoS to another, or 
in simplest form to modify EoS to permit the pres- 
ence of pressure derivatives. 

The simplest way is to introduce the relaxation 
equation for pressure 



tP + P = fip.ait)). 



(4) 



When r = and f{p,a{t)) = P{p) we recover the 
standard EoS. Such an equation may be considered 
as some (dynamical) constraint to usual EoS. Of 
course, the physical sense of such equation (unlike 
to energy conservation law) is not clear at the mo- 
ment although some explanations are given in the 
next section. In daily life, however, there could oc- 
cur similar phenomena where the time change of the 
presure depends on the density. For example, con- 
sider the water. There is a pressure in the steam, 
which is the gas of water. When the density in- 
creases, the molecules of water make drops of water, 
like fog. The pressure of the drops could be ne- 
glected. At high density, the total pressure could 
decrease. The equation seems to express such a 
process. Then if dark energy consists of particles or 
some objects with internal structure, there may oc- 
cur the phase transition like that between steam and 



water drop. At the point of phase transition, since 
the system becomes unstable, the pressure may be 
governed by a equation like 

We prefer to measure the relaxation time t in 
terms of Hubble function H, i.e., consider tH = 
^ ~ const. In this case it is convenient to use a new 
variable x = In terms of x the expression tP 

simplifies as 



A , dP 
tP = ^x— 
ax 



(5) 



and one obtains the pair of relaxation type equations 
for p and P 



^ dp 

rTx^p=-^^ 

dP , , 

^.-+P = /(p,.). 



(6) 
(7) 



Extracting P from the first equation and inserting it 
to the second one we obtain the second order, master 
equation for the energy density p 

"^|S+-$f4+iV7[P + /(P^-)]=0 (8) 



dx 



This is new, dynamical equation to energy-density 
which is compatible with energy conservation law. 

As the explicit example, let the function / be of 
the form (of course, more complicated choices may 
be considered) 



/(P,a;) = -p + 7(x)(p-pc), 



(9) 



where pc = const is some critical value of the energy 
density, and ^{x) = 70 + ax™. When Pc — ^ and 
a = we recover the standard linear EoS 



F=(7o-l)p. 



(10) 



The equation for p can be reduced to the Bessel 
equation and the solution is of the form 

p = p,^x" {CxJ^iax^) + C2 J-.(aa;^)] . (11) 

Here 



1 



A 



a = a 



m/2 
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Using ©, we find 

P = -Pc 

- + 1) x^ [C\JAdx^) + C2 J-.(ax^)] 

[Ci {j,-i{dx^) - j,+i{dx^)) 



(12) 



Xd 



,<7 + X 



+C2 {J-^-i{dx^) - j-^+i{dx^))] 



(13) 
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When X is large, p behaves as an oscillating function 



2v + l 



+C2 cos I ax^ 



-2v+l 



(14) 



Since cr — A/2 = —2 — (3/4)m — 1/^, if we naturally 
assume that m and ^ should be positive, the second 
term damps with oscillation. Then p goes to a con- 
stant p pc- On the other hand, for large x, P 
behaves as 



P 




COS ax 



2v- 1 



+C2 cos ax 



-2v+l 



(15) 



Since a + A/2 = —2 — m/A — 1/^, when m and ^ 
are positive, the second term damps with oscillation 
again and P goes to a constant P ^ —Pc- Then the 
effective EoS parameter w = P/ p goes to —1, which 
corresponds to a cosmological constant. 

The Bessel function is quasi-oscillating and we ob- 
tain an infinite number of epochs, in which p, P, H 
and a are also quasi-oscillating. In other words we 
have an infinite number of points in which the de- 
celeration replaces the acceleration and vice-versa. 

The presence of pc can guarantee that p is positive, 
thus, is also positive. Nevertheless, P can change 
its sign, and this phenomenon can mimic the dark 
energy effect. 

When Pc — a — the equation becomes of the 
Euler type, and the solution is also very simple. 



By using the conservation law (Q, one can rewrite 
(|17() in a form similar to 

tP + P 

= g {p, a) + tH (-2,^^^^ {p + g {p, a)) 



dp 



dg (p, a) 

7) ' 

oa 



(18) 



^H'^p 



(19) 



When other contributions to the energy density can 
be neglected, the first FRW equation looks as 

3 

Then Ea. H18|) can be rewritten as 
rP-l-P 

= g [P, a) + TK^ - 1-3 — — — (p + g [p, a)) 



dg (p, a) 



(20) 



By comparing H2U|) with Q), we may identify 

/ \ , [P f ^dg{p,a) 
= g [P, a) + - 1-3 — — — [p + g [p, a)) 



+ 



dg {p, a) 



(21) 



This shows the relation between standard (gener- 
ally speaking, inhomogeneous EoS,Sj) and relaxation 
equation for pressure. 



IV. GENERALIZED INHOMOGENEOUS 
EQUATION OF STATE 



III. THE RELATION WITH STANDARD 
EQUATION OF STATE. 

We now consider the relation with the standard 
EoS. Let us start with the scale factor dependent 
EoS: 



(16) 



Then we have 
d 



-t/ 

1 

r 
1 



-p + p 



, X , dg (p, q) . dg (p, g) 

(p, a) + ^^p+^^ai7. (17) 



As it was indicated above, there is a possibility 
that the EoS contains P or even higher time deriva- 
tives of pressure. More generally, the EoS could de- 
pend on H or H (inhomogeneous EoS 6]) like 



Uip,P,P,H,H) = 



(22) 



Note that many effective dark energy models like 
brane- worlds, modified gravity and string compacti- 
fications have such a form ( for very recent example 
compatible with observational data, see and ref- 
erences therein). As particular example, one may 
consider 

C/(p, P, P, H, H) = ^+ (f - + ■ 

(23) 
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Here W{p) is a proper function of the energy density 
p. Using the energy conservation law one gets 



P = Wip) 



(24) 



If p is regarded as a coordinate, Ea. (|24|l has a form of 
Newtonian equation of motion of the classical par- 
ticle with the "force" W. For example, if W is a 
constant = , we find p behaves as a coordinate 
of the massive particle in the uniform gravity: 



^0 + ^2 , 



(25) 



Here to and Cq are constants of the integration. As 
an another example, we may consider the case of the 
harmonic oscillator: 

W{p) = -u;^p-po) . (26) 

Then an oscillating energy density follows 0, 0| : 

p ^ po + Asm [tot + a) . (27) 

If other contributions to the energy density may be 
neglected, by using the first FRW equation ifT^ . we 
find the behavior of the Hubble rate, for H25|) . 



(28) 



As ij < when t < t^, and H > when t > to, there 
is a transition from non-phantom era to phantom 
one at t = tf). For H27|) . we have oscillating H: 



H = —yJpo + Asin{ujt + a) . (29) 
v3 

When we neglect the other contributions to the en- 
ergy density and pressure, we also have 



2 . 

—H = p + p 



(30) 



Combining H3U|) with (|19|l . one may define the effec- 
tive EoS parameter W(,g by 



Wcff 



-1 



2H 



Hence, for ^ 



V3K(ifi(t-io)' + co) 



3/2 



(31) 



(32) 



which surely crosses WcS = — 1 when t ^ to. On the 
other hand, for H29|l . one gets 



WcB = -1 - 



2Aijj cos (cut + a) 



^/3K {po + A sin [ut + a)) 



3/2 



, (33) 



which oscillates around WcS = —1 as in |l3l |. 

As an another example, we consider the EoS 

P-3H{p + P) = U{H) . (34) 

Here U (H) is a proper function of the Hubble rate 
H. Then by using (Q, one arrives at 

p + P = U{H). (35) 

In a simplest case, U{H) — 0, it follows 

p + P = c (c : constant) . (36) 

When the other contributions to the energy density 
and pressure are neglected, because of (|5n|) . we find 
H is constant and 



H = -^t 
2 



As an another case, we may consider 
U{H) ^ . 



(37) 



(38) 



Here is a constant. Then combining l|35|l . H30|l . 
and H38|l . we find 



H = -uj^H , 



(39) 



which is the equation typical for the harmonic oscil- 
lator in classical mechanics. Hence, the oscillating 
Hubble rate is obtained 



H = Ho sin {cut + a) 



(40) 



Here Ho and a are constants of the integration. 
Thus, we demonstrated that inhomogeneous gener- 
alized EoS (linear in the pressure derivative) leads to 
the interesting accelerating (often oscillating) late- 
time universe. 



V. THE EQUATION OF STATE 
QUADRATIC ON THE PRESSURE 
DERIVATIVE 

In this section, as an immediate generalization, 
the case that the EoS is not linear on P but 
quadratic is considered. 

Let the equation to pressure and its derivatives 
looks like an energy in the classical mechanics: 



E=-P^ 



V{P) 



(41) 
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Here E is a constant but as it is an analogue of the 
energy, it is denoted as E. We should note that E 
does not correspond to real energy in universe. This 
may be also considered as implicit form of EoS. 
First example is 



V{P) = dP 



(42) 



with a constant d. Then by the analogy with the 
classical mechanics, we find 



1 2 

E = ^vl + dpo . 



(43) 



Here vq and po are constants. 

In case that other contributions to the total energy 
density are large, as in the early universe, the Hubble 
rate H could not be so rapidly changed. Then we 
may assume that the Hubble rate H could be almost 
constant H = Hq. Using one obtains 



P = 



Poe 
d 
~2 



-3Hot 

2 



-Vo 



2t 
ml 
t 



3HoJ 



Po 

3Ho I 3Ho 



(44) 



Here po is a constant. The explicit form of (inhomo- 
geneous) EoS may be found combining two above 
equations. 

On the other hand, we may also consider the case 
that the other contributions to the energy density 
and pressure are neglected as in late-time or future 
unverse. Then deleting p from |^ and (|19|l . we have 



H+^H'+--P = 



(45) 



For igSI), Eq.|g3) admits the solution 

H ^hot + hi , (46) 

when 

d = 3/io I vq = —Shohi , 

Po = ^ho - hi . (47) 

For H46|l . the effective EoS parameter WcS defined by 
(pTT|l has the following form: 



WcS = -1 



2hQ 



3{hot + hiY 



(48) 



which goes to —1 when t goes to infinity. Hence, the 
emerging universe seems to be the asymptotically de 
Sitter one. 

Second example is 



Then we have 



V{P) = ^co^P' 



P = Aam{ujt + a) 
E = A^u? . 



(49) 



(50) 



where A and a are constants. Then in the case 
that other contribution to the total energy density 
is large, as in the early universe, the Hubble rate H 
could be almost constant H = Hq, wc find 

-3^0* 



P = Poe 



'UJ cos [ujt + a)) 



(51) 



with a constant po- This corresponds to de Sitter 
universe. 

On the other hand, when other contributions to 
the total energy density can be neglected, as in the 
late-time universe, by using (|45ll . one gets 



dPa'/^ 3 
df^ ^ 4 



+ -K'^Asin{u}t + a)a^/^^Q. (52) 



By defining a new variable s 

s = Lot + a + — , 
one obtains a kind of Mathieu equation: 



whose solution is given by 



cos s a 



3/2 



(53) 



(54) 



^3/2 



a"' ' = Cn cos{nt) + s„ sm{nt) . (55) 



Here the coefficients c„ and s„ are given by recur- 
sively solving the following equations: 

Co = c , ci = 

-n'^Cn + I (c„_i + c„+i) = (?i > 1) , 
2 

Si = S , S2 = --S , 

-n^Sn + I {Sn~l - Sn+i) =0 (n > 2) , 



3k^ 
4^ 



(56) 
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Hence, a has a periodicity 1/oj. In the expression 
()55|l . a is not always positive. Then physically the re- 
gions where a'^/^ is not negative could be allowed and 
the points a = could correspond to Big Bang/Big 
Crunch/Big Rip^. 

We should note that the expressions of po in 
and H51|) are not always positive. Then only the 
period(s) where po is positive could be allowed in 
the real universe. 



VI. COUPLING WITH THE MATTER 



which gives, as well-known, 



H 



3(l + u)„) 
t 



(63) 



On the other hand if ct — Ai^ < —3(1 + Wm) < 0, 
Hubble rate is 




which gives 



1/2 



,(<T-A^)/2 



(64) 



No direct interaction between dark energy 
and matter 



Let us now include the matter. For simplicity, we 
consider the matter with constant EoS parameter 
Wm so that the matter energy density pm is given by 



P = Po 



a{to) 



(57) 



In case of the total energy density is given by 

Ptot = Pc+x" [CiMax^) + C2J^Aax^)]+Pox'''^^^'' 

(58) 

and the Hubble rate H is given by 

-t-pox-3(i+-™))}'^' . (59) 



In future, x becomes large, then the Hubble rate H 
goes to a constant (with oscillations): 

(60) 

which tells WeS —1. On the other hand, in the 
early universe, x should be small. Hence, one finds 



1/2 



(61) 



If — 3(l + m,Ti) < (T — A, the contribution from matter 
becomes dominant and Hubble rate is 



H = 



4f} 



1/2 



-3(l+tu„)/2 



(62) 



H 



(7 — \u 



(65) 



By comparing H65() with (|63|l or it follows that 
the effective EoS parameter is given by 



a — Xv 



(66) 



For the model in solving (|^ . we find the 

^-dependence of p. Then the FRW equation gives 



= Pit) 



Po 



a (to) 



-3(l+tu™) 



(67) 



For the case (|25|) . when t is large enough, the second 
term in the r.h.s. of H67|) could be neglected and we 
will obtain H28|l . If cq in H25() is small enough, when 
t to, the second term in (|67|l could be dominant 
and we may obtain H63|) . For the case H29|) . in the 
early universe, where a is small, the second term in 
H67() could be dominant and one obtains (|63|l . again. 
Especially for the dust w„i = 0, we find H ~ 
that is, a ^ ti . In the late time universe, the first 
term could be dominant and one gets (j^ . 

Three years WMAP data are recently analyzed 
in Ref. 16], which shows that the combined analysis 
of WMAP with supernova Legacy survey (SNLS) 
constrains the dark energy equation of state wde 
pushing it towards the cosmological constant. The 
marginalized best fit values of the equation of state 
parameter at 68% confidance level are given by 
— 1.14 < Wde < —0.93. In case of a prior that 
universe is flat, the combined data gives —1.06 < 
Wde < -0.90. 

In our models, as shown in p6|l . H17|l . 148|) . and 
(|^ . the effective EoS parameter is WeS ~ — 1 and 
there is no contradiction with the above WMAP 
data. We should also note that when mater is cou- 
pled, we find WcS ^ Wm in the early universe, as in 
ijB^ . Thus when Wm < —1/3, there should occur 
the transition from deceleration to acceleration. 
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B. Dark energy interacting with matter 

Generally speaking, the matter interacts with the 
dark energy. In such a case, the total energy den- 
sity ptot consists of the contributions from the dark 
energy and the matter: ptot = P + Pm- If we define, 
however, the matter energy density pm properly, we 
can also define the matter pressure p„i and the dark 
energy pressure p by 



P 

3H 



P = Ptnt ^ Prr 



(68) 



Here Ptot is the total pressure. Hence, the matter 
and dark energy satisfy the energy conservation laws 
separately, 



p,i 



Pm) = 



p + 3H{p + P) = . 

(69) 

In case, however, that the EoS parameter w„i for 
the matter is almost constant, one may write the 
conservation law as 

Pjn+iH{l+Wm)Pm^Q , (70) 

and therefore for the dark energy 

p + 3Hip + P) = -Q , (71) 

so that the total energy density and the pressure 
satisfy the conservation law: 



Ptot + 3i7 (ptot + Ptot) Pm^O 



(72) 



In H7()|l , Q expresses the shift from the constant EoS 
parameter case. 

As an example, we consider the case that Q is 
given by a function q ^ q{a) as 

Q = Haq'{a)pm ■ (73) 
Combining (|7^ with (|7n|) . one gets 

P,n = p,„oa-'^'+'""^e«('') . (74) 

Here pmo is a constant of the integration. Hence, the 
conservation law is modified, through ^71^ as 



p + 3Hip + P)^ -p,„oi?a-^'+"""'9'(a)e''('^' 
and © is also modified as 

L^ + p = -P-s{x) . 

6 ax 

Here 



(75) 
(76) 



S{x) ^ {a{to)x)-^'+''"-^ q' {a{to)x) e^f^f*"'-) . 

(77) 



Note that Eq.® is also modified: it now contains 
the inhomogeneous terms: 

2 (Pp , ^dp f , 1 \ 3 , 



3a; 



dS{x) 
dx 



-S[x) 



dx 



(x^'^S{x) 



(78) 



Let a (special) solution oi ^ he p = Ps{x). Then 
in case of (O with 7(0;) = 70 + ctx™, the general 
solution corresponding to (|ll|l is given by 

p = p,{x) + [C^J,{ax^) + C2 J-,(ax^)] . (79) 

where Pc should be included in Ps{x). It is also noted 
that the initial conditions are relevant to determine 
Ci and C2 but irrelevant for Ps{x). As an example, 
we can find 

p,(x) =pc + poa;'' , (80) 
with constants po and 77 when 

p9(o) ^ p<30 

(772 - 3r;/4 + r;/e + 370/O ''a''+3(i+-") 



e^^ ' = e 
Pa 

PmO 



(77+1/0 iv + Sil + Wra)) 

(3a/Oao""'a'"+''+3(i+''"-) 

" (771 + 77 + 1/0 (771 + 77 + 3(1 + W,n)) 

3.soay«a-i/«+3(i+-"') 



(-1/^ + 3(1 + ^;™)) ■ 
which gives 

' (772 - 377/4 +77/e+ 370 /e) 



(81) 



S 



Po 
3 



3ax™+'' 



(m + 77 + 1/C) ^ 



SqX 



(82) 



In H81|) and (|82|l . go, so are constants and oq = a(io)- 
In case that 

77+3(l+w„0, 77i+?7+3(l+w,„), -l/^+3(l+u;,„) < , 

(83) 

we find e''^"-' e'^° when a becomes large, that 
is, in the late time universe. Thus, pm — > 
Pmoe''"a~(2+'^'"'"^. Furthermore if 77 < 0, we find 
Ps — * Po, that is, H goes to a constant, which may 
lead to the asymptotically deSitter space. Clearly, 
for more complicated coupling Q, more sophisticated 
accelerating cosmology may be constructed. 
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VII. DISCUSSION 

In summary, we discussed the constrained EoS 
for cosmic fluid where the relaxation equation for 
pressure is introduced. It is shown that such 
EoS is equivalent to usual inhomogeneous EoS Q 
which contains scale factor dependent terms. Sub- 
sequently, the generalized inhomogeneous EoS with 
time derivatives of pressure is presented. For the 
number of explicit examples, the accelerating dark 
energy cosmology as follows from such EoS cosmic 
fluid is constructed. It turns out to be the asymptot- 
ically de Sitter universe or oscillating universe with 
long accelerating phase and transtion from decelera- 
tion to acceleration. The consistent coupling of such 
constrained EoS dark fluid with matter is discussed. 
It is shown that emerging FRW cosmology may be 
consistent with three years WMAP data. 

Of course, there are many ways to generalize the 
EoS for cosmic fluid and to investigate the corre- 
sponding impact of such generalization to dark cos- 
mos. The physics behind such generalization re- 
mains to be quite obscure (as dark energy itself 
and its sudden appearence). At best, this may be 



considered as some phenomenological approxima- 
tion. Nevertheless, having in mind, that most of 
modern attempts to understand dark energy includ- 
ing strings/M-theory, brane- worlds, modifled grav- 
ity, etc lead to effective description in terms of cos- 
mic fluid with unusual form of EoS, it turns out to 
be extremely powerful approach. From another side, 
the reconstruction of the cosmic fluid EoS may be 
done for any given cosmology compatible with obser- 
vational data which may flnally select the true dark 
energy theory. 
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